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A theoretical description of radiation-matter coupling for semiconductor-based photonic crystal 
slabs is presented, in which quantum wells are embedded within the waveguide core layer. A full 
quantum theory is developed, by quantizing both the electromagnetic field with a spatial modulation 
of the refractive index and the exciton center of mass field in a periodic piecewise constant poten- 
tial. The second-quantized hamiltonian of the interacting system is diagonalized with a generalized 
Hopfield method, thus yielding the complex dispersion of mixed exciton-photon modes including 
losses. The occurrence of both weak and strong coupling regimes is studied, and it is concluded that 
the new eigenstates of the system are described by quasi-particles called photonic crystal polam- 
tons, which can occur in two situations: (i) below the light line, when a resonance between exciton 
and non-radiative photon levels occurs (guided polaritons), (ii) above the light line, provided the 
exciton-photon coupling is larger than the intrinsic radiative damping of the resonant photonic 
mode (radiative polaritons). For a square lattice of air holes, it is found that the energy minimum 
of the lower polariton branch can occur around normal incidence. The latter result has potential 
implications for the realization of polariton parametric interactions in photonic crystal slabs. 

PACS numbers: 71.36.+C, 42.50.Ct, 42.70.Qs, 71.35.-y 
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I. INTRODUCTION 



The capability of engineering electron and photon 
states through spatial confinement leads to a control of 
radiation-matter interaction and to a number of interest- 
ing results in the field of solid-state cavity quantum elec- 
trodynamics (QED)ii For example, the use of quantum- 
well (QW) excitons embedded in high-finesse semicon- 
ductor microcavities (MC) of the Fabry-Perot type has 
allowed to observe a modification of spontaneous emis- 
sion (weak coupling regime )2 as well as the occurrence of 
a vacuum Rabi splitting (strong coupling regime) 4ii The 
latter effect arises when the radiation-matter coupling 
energy overcomes the damping rates of QW excitons 
and MC photons. Under such conditions, the elemen- 
tary excitations of the system should not be described 
as barely excitonic or photonic, but rather as mixed 
radiation-matter states, called MC polaritons~&Z&& The 
physics of exciton-polaritons in bulk materials had been 
introduced in the fiftie s 10 ' 11 and is now a textbook 



topic 
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The confinement of both excitons and pho- 



tons with MC-cmbcddcd QWs leads to polariton states 
that are more easily observed, even at room temperature, 
and in fact MC polaritons have been a major research 
topic in the nineties. In 2000, the demonstration of co- 
herent parametric processes due to polariton-polariton 
scatterin g 14 ! 15 opened a new area of research in the 
physics of MC polaritons and nonlinear optics. More re- 
cently, the experimental demonstration of Bosc-Einstein 
condensation of MC polaritons with II- VI materials^ 
brings a new exciting development and demonstrates the 



potentialities of radiation-matter interaction in the pres- 
ence of electron and photon confinement. 

Until now, little has been done to study radiation- 
matter interaction when high-quality QWs are combined 
with new routes for the tailoring of electromagnetic field 
in photonic crystals (PhC). These structures, which are 
characterized by a periodic modulation of the dielectric 
constant along one, two or three dimensions (ID, 2D, 
3D)jiZi allow for an unprecedented control over the prop- 
agation and confinement of light at optical wavelengths. 
In particular, two-dimensional PhC embedded in a pla- 
nar dielectric waveguide, commonly known as ■photonic 
crystal slabs, are receiving much attention because they 
allow for a 3D control of light and are easily realized 
at sub-micron length scales . 18 d 9 i 20 i 21 In these systems, 
light propagation or confinement is controlled by the PhC 
structure in the 2D plane and by the dielectric disconti- 
nuity provided by the slab waveguide in the vertical direc- 
tion. When QWs are embedded in a PhC slab, additional 
flexibility and new effects related to exciton-photon cou- 
pling are to be expected. 

In this work we analyze radiation-matter interaction 
in PhC slabs, by considering the effects of coupling be- 
tween QW excitons and photonic modes. A main goal 
is to study the conditions leading to the strong cou- 
pling regime and to point out the possibilities offered 
by PhC for strong coupling engineering. In the past, 
both experimental^ 2 - and theoretical^ investigations of 
strong coupling regime in organic-based ID PhC slabs 
have been reported, the theory being based on a semi- 
classical treatment. The same groups also reported a 



joint work on 2D PhC slabs with organic active medium, 
in which the theoretical treatment was again based on 
a semiclassical solution of Maxwell's equations within a 
scattering matrix formalism* 2 ^ Exciton-polaritons in ID 
PhCs (ideal Bragg mirrors) in which one of the con- 
stituent media has a strong excitonic character were the- 
oretically studied by using a fully classical approach, 25 
Other theoretical works focused on polaritonic gaps in 
so-called phonon-polaritonic PhCs^ or on the proper- 
ties of 3D arrays of quasi-zero-dimensional excitons (e.g., 
quantum dots) coupled to Bloch modes of the same peri- 
odic array of dielectric objects^I More recently, a purely 
2D square lattice array of rods in a dielectric host medium 
has been theoretically studied, in which one of the two 
constituent media has a frequency-dependent permittiv- 
ity with a pole at a resonant frequency^ Combining 
the effects of exciton confinement in QW structures with 
the 2D periodicity of the surrounding dielectric environ- 
ment has been only studied in relation to the influence 
of Coulomb interaction on the longitudinal part of the 
electromagnetic field for deeply patterned PhCs with an 
unpatterncd QW underneath the PhC region^ No ex- 
perimental or theoretical work has been reported so far, 
to the best of our knowledge, on excitons or exciton- 
polaritons in semiconductor-based PhC structures that 
are the subject of this work. 

In the present paper, the theoretical problem of 
exciton-photon coupling in PhC slabs is tackled by 
using a fully quantum-mechanical formalism for both 
photons and exciton states, which is described in de- 
tail. With respect to previous works on the quan- 
tum theory of exciton-polaritons in semiconductor nano- 
structureS ) 30 ' 31 i 32 ' 33 i 34 ' 35 the non-trivial spatial depen- 
dence of the dielectric constant is taken into account 
by quantization of the electromagnetic field in a non- 
homogeneous mcdium . 36 ' 37 i 38 ' 39 The total hamiltonian 
for the coupled exciton-photon states is derived and di- 
agonalized numerically to obtain the eigenenergies of the 
mixed modes. It is shown that polaritonic effects are 
typically stronger in PhC slabs than in MCs, due to the 
better field confinement provided by total internal re- 
flection in the slab. As a consequence, larger vacuum 
Rabi spitting is found at exciton-photon resonance in 
strong coupling regime, which occurs when the exciton- 
photon coupling energy is larger than the intrinsic radia- 
tive linewidth of the photonic mode and the QW exciton. 
In this scenario, new quasi-particles form in the PhC slab, 
which we call photonic crystal polaritons. Preliminary re- 
sults of the formalism presented here have been published 
as conference proceedings ; 40 ' 41 ^ 42 and limited to ID or 
triangular PhC lattices. In addition to the full theoret- 
ical formulation, we present here new results for a PhC 
slab made of a square lattice of air holes in a high-index 
suspended membrane. For a specific PhC slab design, a 
polariton angular dispersion is found which has a mini- 
mum close to normal incidence with respect to the slab 
surface. Such result paves the way for new exciting de- 
velopments in the investigation of hybrid semiconductor 
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Figure 1: (color online) Schematic view of a two-dimensional 
photonic crystal slab of thickness d and lattice constant a, 
with a single quantum well grown within the core layer. 



structures exhibiting both photonic and electronic band 
gap characteristics, as it will be discussed. 

The paper is organized as follows. In Sec. |TT] we 
schematically describe the system under investigation. 
In Sec. IIIII we give a detailed account of the second- 
quantization procedure leading to the complete exciton- 
photon hamiltonian in the linear (low excitation density) 
regime and of the diagonalization technique. In Sec. IIVI 
we provide systematic results in both weak and strong 
coupling regimes, and describe the formation of guided 
and radiative polaritons. The relevance of the results in 
the context of parametric processes involving polariton- 
polariton scattering is also discussed. Finally, in Sec. [V] 
we summarize the conclusions of the work. 



II. THE SYSTEM 

The schematic of the model system under investiga- 
tion and the choice of the coordinate axes are given in 
Fig. [T] We consider a high- index planar dielectric waveg- 
uide with semi-infinite cladding layers (typically the up- 
per cladding is air). We assume that a QW of thickness 
Lqw is grown in the core of the waveguide at a cer- 
tain vertical position zqw, measured from the core/lower 
cladding interface. The dielectric material constituting 
the core of the planar waveguide acts also as a barrier 
material for the carriers (electrons and holes) confined in 
the thin QW layer. Typical QW thickness is of the order 
of Lqw — 10 nm, while the dielectric core layer is usually 
between d = 100 and 200 nm thick. The QW exciton is 
characterized by its transition energy, E ex , and by the 
oscillator strength per unit area, f/S£l The effect of 
multiple quantum wells (MQWs) grown in the core layer 
can be taken into account by a total oscillator strength 
given by the sum of single QW oscillator strengths. This 
approach is valid as long as Lqw *C d, and the electric 
field of the relevant photonic mode can be considered as 
uniform along the MQW thickness. We disregard the 
difference in dielectric constants between QW and bar- 
rier media, and assume the core layer to be described by 
a relative dielectric permeability edici- The core layer is 
patterned with a given photonic lattice, e.g. with lattice 
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constant a, down to the interface with the lower cladding 
medium. Thus, exciton center-of-mass eigenfunctions are 
not free as in usual QWs, but are subject to a further 
confining potential provided by the etched air regions. 
As a general remark, it is clear that exciton and photon 
wavefunctions are confined in the vertical direction and 
are subject to effective potentials having the same spa- 
tial periodicity in the 2D plane. Roughly speaking, we 
could say that confined excitons and photons display the 
same dimensionality in this problem, thereby satisfying 
the condition for possible occurrence of quasi-stationary, 
strongly coupled polariton states 2^ 

III. QUANTUM THEORY OF 
EXCITON-PHOTON COUPLING 

In this Section, we detail the full quantum theory of 
radiation-matter interaction for QWs embedded in a PhC 
slab. We start by the second quantization of the electro- 
magnetic field in a non-uniform dielectric environment, to 
which end we follow the treatment given in Ref . l37l. Then, 
we quantize the exciton field by solving the Schrodingcr 
equation for the center-of-mass QW exciton wavefunc- 
tion in a periodic, piecewise constant potential. The so- 
lutions for the non-interacting fields are used to derive 
the coupling hamiltonian. The total hamiltonian is then 
diagonalized, leading to PhC polariton eigenmodes. 

A. Quantized photon field in non-uniform 
dielectric medium: normal mode expansion 

The canonical quantization of the electromagnetic field 
in a non-uniform dielectric medium with relative perme- 



ability e(r) (we assume a unit magnetic permeability) is 
carried out through quantization of the vector potential 
A, defined from the usual relations to the fields 

■<'■«> ->■'). ™ 

B(r,i) = VxA(r,t), (2) 
and which must satisfy the equation of motion 



V x V x A = — e(r)A. (3) 

cr 

We are considering only the retarded electromagnetic 
field (no scalar potential, <t> = 0) and the vector potential 
can be chosen to satisfy the generalized Coulomb gauged 

V-(e(r)A(r,t))=0. (4) 

The second-quantized hamiltonian of the free photon 
field is straightforwardly obtained from the classical ex- 
pression of the total electromagnetic energy 

H e . m . = ^ J [e(r)E(r) 2 + B(r) 2 ] dr , (5) 

which can be derived from the Proca Lagrangian 
density^ and can be expressed in terms of A(r, t) and its 
conjugate momentum II(r, t) = e(r)A(r, t)/c 2 . The field 
operator is expanded in normal modes as 



A(r,t) 



k . n 



1 /2 



kn^-k 



n( r )e 



(6) 



in which a^ n (aun) are creation (destruction) operators 
of field quanta with eigenenergies <x> k „, and the index n is 
a generic band number labelling the corresponding pho- 
tonic eigenmode at in-plane Bloch vector k. In order to 
satisfy Bosc commutation relations for the operators, 

[aun, Ok'n'] = <5k,k'<5ri.n' , [akra, Ok'n'] = [Ofcri' <*Vn'1 = ®> 

(7) 

the correct orthonormality condition for the classical 
functions Akn(r) satisfying Eq. ((3]) ig2L2LSi 

/ e(r)A^(r).A k ^(r)dr = ^<5k,k'W, (8) 
Jv w k„ 



where V is the same quantization volume as in ([5]). This 
implies the following conditions for the electric and mag- 
netic fields: 

/ e(r)E k „(r)-E k ^(r)dr = «5 k , k ,«5 nX , (9) 
Jv 

f B kn (r) • B k <„- (r)dr = Sk,k'<W ■ (10) 
Jv 

The usual Maxwell equations for E and B can be de- 
rived from the equations of motion for the conjugate vari- 
ables A(r) and II(r). From Eq. ((5|) and using the field 
expansion ^ . we finally get the second-quantized hamil- 
tonian of the free photon field in the usual form as a sum 
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of harmonic degrees of freedom. 



-^p' 1 = E ^ Wk ™ ( "kn" k 



(11) 



In order to implement the formalism, we shall need 
specific forms for the classical fields functions. The exact 
solution of Maxwell equations in a PhC slab is a com- 
plicated task, especially concerning quasi-guided modes 
that lie above the light line i 18 ' 19 A recently developed 
guided-mode expansion (GME) approach^ allows us to 
find a convenient solution to this problem after expansion 
of the classical fields on the basis of guided modes of an ef- 
fective planar waveguide. We start from the second-order 
equation for the magnetic field B = H in a source-free 
dielectric medium and for harmonic time dependence, 



V x 



1 

W) 



V x H 



4« 

C z 



(12) 



Due to the in-plane translational invariance implying 
Bloch-Floquet theorem, the magnetic field can be ex- 
panded on a basis in which planar and vertical coordi- 
nates are factorized 



H k (r) = ^^c k , Q (G)h k+G . Q (z)e 



j(k+G)-p 



(13) 



G 



where r = (p, z), k is the in-plane Bloch vector in the 
first Brillouin zone (BZ), G are reciprocal lattice vec- 
tors, and the functions hk+G,ct(z) (a = 1,2,...) are the 
(discrete) guided modes of the effective planar waveguide 
with an average dielectric constant in each layer, calcu- 
lated from the air fraction of the given photonic lattice. 
Thus, Eq. (|12p is reduced to a linear eigenvalue problem 



the second-quantized formulation. Finally, we point out 
that losses of a quasi-guided photonic eigenmode can be 
introduced in the present treatment as an imaginary part 
7kn of mode eigenenergies. Such imaginary part arises for 
modes lying above the light line from out-of-planc diffrac- 
tion losses, which are calculated in perturbation theory 
as detailed in Ref. l4rl 
B. Quantized exciton field in a periodic piecewise 
constant potential 

Exciton confinement in QWs has been widely investi- 
gated in the past (for a review see, e.g., Ref. [431 ). Also, 
second-quantization of the QW exciton field has been ex- 
tensively treated in the literature. Two recent approaches 
leading to an effective hamiltonian of interacting bosons 
starting from the full crystal hamiltonian are given in 
Refs. I47ll48l When the exciton density n ox is low, i.e., 
much smaller than a saturation density n sat ~ 1/(2^0^) 
(where a-io is the 2D exciton Bohr radius) ) 49 i 50 QW ex- 
citons behave as a gas of non-interacting bosons to a very 
good approximation. 

The envelope function F(r c ,rh) of QW excitons can 
be factorized in a wavefunction depending on a ccnter- 
of-mass coordinate R|| and an electron-hole wavefunction 
depending on the in-plane relative coordinate p as 

F(r 0> r h )=F K (R||)/(p,« ej ^), (15) 

where Fk(Rii) = e lK R » for free QW excitons. Consider- 
ing the system of Fig. [TJ the QW layer is patterned with a 
photonic lattice of air holes and the center-of-mass wave- 
function Fk(R||) for exciton motion in the 2D plane is 
not a simple plane wave anymore ^ This wavefunction 
obeys a Schrodinger equation 



E E ^+G,k+G'<W (k + G') = ^c Q (k + G) , (14) 



G' 



2M C > 



mil) 



Fk(R||) = Ek Fk(R||) . 



(16) 



which is solved by diagonalizing the explicit expression 
for the matrix M. 4^ The properties of the specific pho- 
tonic lattice enter in the matrix M. as a Fourier trans- 
form of the inverse dielectric constant in each layer, 
r](G,G') = e _1 (G,G'), the matrix inversion being per- 
formed numerically. After diagonalization, the result- 
ing photonic modes can be classified according to their 
band index, n, and their in-plane Bloch vector k. It 
should be noted that the classical fields H k „(r) and 
E k „(r) = ic/(we(r))V x H k „(r) calculated by the GME 
approach automatically satisfy the orthonormality condi- 
tions ([9l fl0| , thus they constitute a very convenient set for 



where M cx = m* + is the total exciton mass and the 
potential ^(R||) = in the non-patterned regions, while 
y(R||) takes a large value in the air holes. Equa- 
tion (|16p can be solved by plane-wave expansion, 



F K (Rt|)=E^( K 



G 



G)e i(K+G).R„ ^ (17) 



where the exciton wavevector K is now restricted to the 
first Brillouin zone and G are the same reciprocal vectors 
of the photonic lattice. The resulting equation in Fourier 
space is 



J 



E 

G' 



r ft 2 |k + Gj 

2M ex 



#G,G' + V(G — G') 



F{K + G') = E K F(K + G) 



(18) 



5 



where the Fourier matrix V(G, G') = V(G — G') has 
expressions similar to those for e(G, G') in the photonic 
problem. Equation (|18| is solved numerically, yielding 
quantized center-of-mass levels in the periodic potential. 
By this procedure, the exciton levels are labelled by the 
same quantum numbers as the electromagnetic modes, 
namely a Bloch vector K and a discrete index v. The 
exciton energy is written as E^} = E cx + Ejc u , where 
E ox is the bare QW exciton energy (which is given for the 
specific QW parameters) , while E\^ v is the center-of-mass 
quantization energy in the in-plane potential V(R||). It 
turns out that E^ ~ 1CP 2 — 10 _1 mcV for typical 
PhC slab patterns, depending on a and r/a. Now we 
may introduce exciton creation (destruction) operators 
for center-of-mass eigenmodes, b^ v (&kj/)- In second 
quantization, the hamiltonian of the bare exciton field 
is finally given by 

^cx^E^K^L^K,, (19) 
K,i/ 

where for n ex 1 /a\ D the excitonic operators obey com- 
mutation relations 

(20) 



where thq is the free electron mass, Tj (jpj) are the posi- 
tion (momentum) variables of the QW electrons, and the 
sum runs over all the electrons in the system. In Eq. (|21[) . 
we have retained both A • p and p • A terms, because the 
generalized Coulomb gauge ^ does not allow in gen- 
eral for the commutation of these two operators. Taking 
into account the presence of a non-local potential, we 
can write this hamiltonian in second-quantized form as 
Hi = + £Tj 2 , where the two contributions ar o 52 ' 53 



N 

= ^E A fo)-fc>^ex], (22) 
2 N 

^ = -^D^)-^^)-^-], (23) 



and the operator Vj = Pj/m = (ifr)^ 1 [Yj, H ex ] directly 
follows from Heisenbcrg equation of motion for ij . Intro- 
ducing exciton operators and resolving the commutator, 



Broadening of the exciton spectral lines can be taken 
into account phenomenologically by introducing an imag- 
inary part of mode energies 7 0X , taken as a K- and v- 
independent parameter. 



C. Exciton-photon interaction hamiltonian 



The interaction between exciton and photon states 
must conserve the Bloch vector, i.e., K = k (mod- 
ulo a reciprocal lattice vector). Thus, we will use the 
same notation for exciton and photon wave vectors. 
However, in general a photonic mode with band index 
n couples to exciton center-of-mass levels with any v. 
The interaction is determined by a matrix element of 
the full hamiltonian, as first shown in Refs. fl^flll for 
bulk exciton-polaritons and later extended to quantum- 
confined systems . 30 ' 33 ! 34 ' 35 The classical minimal cou- 
pling hamiltonian is given by 



(21) 

I 

one can write, e.g., for the first term 

= E ^ I E Mr 3 ) ■ r 3 10) bl + h.c. , 

k,i/ j 

(24) 

where is the many-body exciton wavefunction while 
0} is the crystal ground state. Expanding the vector po- 
tential as in Eq. ([6]), and expressing also Hj in terms of 
exciton operators, after some manipulation the two terms 
of the second-quantized interaction hamiltonian are ob- 
tained in the form 

= i ]T C w (a k „ + al kn )(l^ - b{ v ), (25) 

k,n,i/ 

A <2) = E D *vnn' (6-kn + aLjOW + oL^Jftfi) 
k,z/,n,n' 

The coupling matrix element C kn „ between exciton and 
photon states at a given k is calculated as^ 4 - 

1 /2 

C^nv = I ^ j (* kt/ | 2^ A k „(r,) • r 3 10) , 

(27) 



N 



N 



Hl = ~ W ^ ' ^ + Vj ' A(r ' ? ' )} 



J'=l 



2moC" 



Ei A ^-)i 



i=i 
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while D klW = C^C^JE^ 

The integral in Eq. (j2"T)) can be expressed in terms the 
oscillator strength / of the excitonic transition, which is 
generally defined as^ 3 - 



(ex) 



,-|o)| 2 



(28) 



where Q cx = E^y 1 / h and e is the polarization unit vector 
of the exciton. For a QW exciton the oscillator strength 
per unit area is calculated as 



h 
s 



2ninQ 



£ cx i 



e • r c 



f(p = 0,z,z)dz 



(29) 



in which r cv = (w c o|r|u v o) is the dipolc matrix clement 
between the single-particle Bloch functions in the valence 
and conduction bands of the bulk crystal. We are consid- 
ering the ground-state heavy-hole (HH) exciton, whose 
optically active states are polarized in the xy plane and 
are doubly degenerate. Thus, HH QW excitons prefer- 
entially couple to TE-like modes in the photonic struc- 
ture. Expressing the exciton wave function in terms of 
the QW envelope function (fT~5|) . it is easy to see that the 
coupling energy (|27p depends on the oscillator strength 
per unit area, Eq. ([29]) . as well as on the spatial overlap 
between the exciton center-of-mass wavefunction and the 
photonic mode profile in the QW plane. Assuming the 
mode profile to be uniform along the QW thickness, we 
can express the coupling matrix element as 5 ^ 



C 



nh 2 e 2 f 
toq S 



1/2 



J e-E kn (p, zq W )F^(p) dp, 

(30) 



where we approximated E^^ ~ hw^n, assuming close- 
to-resonance coupling^ We notice that for MQWs the 
coupling energy ([30)) can be multiplied by y/Nqw, where 
Nqw is the effective number of QWs coupled to the pho- 
tonic mode. Within our GME formalism, the integral in 
Eq. ([30]) can be straightforwardly calculated in Fourier 
space and we finally get 



toq S 



1/2 



G 



e • Ek+G,ra(2Qw)-Fk+G,fc 



(31) 

where Ei c +G,n( 2; Qw) is the Fourier transform of the mode 
electric field at the QW vertical position^ Thus, all pa- 
rameters of the interaction hamiltonian are obtained in 
terms of the electric field coefficients calculated by the 
GME method and of the exciton coefficients correspond- 
ing to quantized ccntcr-of-mass levels. 



D. Diagonalization of the total hamiltonian 

The full quantum hamiltonian describing the coupled 
QW exciton and PhC slab modes is given by 



P h 



(32) 



where i? p h, H cx and the two interaction hamiltonians 
are explicitly given in Eqs. {11]), {!§]), and ([2"5|) - ([2"o]) . 
In the bare photonic dispersion, we do not consider 
the zero-point energy term. Not surprisingly, Eq. (|32p 
has a formal analogy with the second-quantized exciton- 
photon hamiltonian in bulk,-? planar MCs^ or pillar 
MCsi^ 5 . Such hamiltonian is valid in the linear regime, 
i.e., under low excitation. Non-linear terms includ- 
ing exciton-exciton scattering and saturation of exciton- 
photon couplin g 49 ! 50 are not considered here. 

The total hamiltonian of the exciton-photon coupled 
system is diagonalized by using a generalized Hopfield 
transformationi 10 i 30 ' 34 i 35 New destruction (creation) op- 
erators ?3 k {pi) are defined as 

Pk = W kli flkn+^] Zk^k^+J^ yknalk ra +yi z k^l k „ , 

nun V 

(33) 

which still satisfy Bosc commutation relations 

\p*,pl]=h,v, [Pk,Pk'] = [P k ,pl]=0. (34) 

The condition for the total hamiltonian to be diagonal in 
terms of p k , p^, is 



[Pk, H u 



KlkPk ■ 



(35) 



The transformation, which leads to a non-hermitian 
eigenvalue problem, applies also to a hamiltonian that 
includes dissipative terms. This is the present case when 
the imaginary part of the frequency for quasi-guided pho- 
tonic modes, as well as the exciton lincwidth arising from 
non-radiative processes, are included in the terms -ff p h, 
H cx of Eq. ([3"2")l . The eigenvalue problem can be written 
in the form 



(36) 



where M k is the generalized Hopfield matrix at a spe- 
cific wave vector k, and v k is a generalized vector con- 
stituted by the expansion coefficients of Eq. ([55]) (see 
Appendix). Diagonalization of Eq. (|36[) gives directly 
the complex eigenenergies ?if2 k corresponding to mixed 
excitations of radiation and matter. Depending on the 
interplay between exciton-photon coupling and their re- 
spective losses, the system can be either in a weak or 
in a strong-coupling regime. In the latter case, /if2 k cal- 
culated for any k in the first BZ gives rise to the full 
spectrum of photonic crystal polaritons. 



IV. NUMERICAL RESULTS: THE SQUARE 
LATTICE PHC SLAB 



In this Section we apply the present theory to a square 
lattice of air holes in a symmetric PhC slab with air 
claddings. We consider a GaAs membrane with a dielec- 
tric constant edici = 12.46, typical of this material in the 
near infrared. A schematic representation of the direct 
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Figure 2: (color online) (a) Square lattice of air holes: direct 
and reciprocal lattices with Brillouin zone and main symmetry 
points, (b) photonic mode dispersion (even modes, <j xy = +1) 
when the lattice is embedded in a high-index (ndid = 3.43) 
photonic crystal membrane. Parameters are: r/a = 0.34, 
d/a = 0.3. The first few modes are labelled by a band number. 
Dashed lines represent the light dispersion in air. 



and reciprocal lattices is shown in Fig. HJa). The pho- 
tonic lattice is characterized by the inter-hole separation, 
a, and the hole radius, r. In reciprocal space, the BZ is 
evidenced and the main symmetry points T, X, and M 
are defined. In Fig.[^b) we show a typical photonic band 
diagram calculated by using the GME method, 46 for a 
PhC membrane with thickness d/a = 0.3 and hole radius 
r/a = 0.34. The dispersion is shown in dimcnsionlcss 
frequency units [uja/{2irc) = a/A], and only for modes 
with even parity io xy = +1) with respect to a horizontal 
mirror plane at the middle of the slab. Such modes have 
mainly in-plane polarized electric field (also defined TE- 
like modes) ) 19 i 46 and they are dominantly coupled to HH 
exciton states in the QW. No higher-order modes of the 
slab waveguide are present for this structure in the fre- 
quency range of Fig.[^b). The first five photonic bands 
arc labelled with integer numbers. Mode 1 and part of 
mode 2 lie below the air light line and are therefore truly 
guided. In principle, such modes have infinite lifetime, 
and can leak radiation out of the slab plane only due 
to disorder in the photonic structure ) 57 ' 58 which we ne- 
glect in the present paper. On the contrary, the modes 
are quasi-guided when their dispersion falls above the 
light line in the first BZ. These modes are coupled to 
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Figure 3: (color online) Guided polariton dispersion for a 
photonic mode interacting with a QW exciton at E cx = 1.485 
eV, close to the BZ edge (X-point); the uncoupled mode dis- 
persions are shown with dashed and dot-dashed lines, respec- 
tively. Parameters of the photonic structure are: r/a = 0.34, 
d/a = 0.3. Panel (a): coupling to the first photonic mode, lat- 
tice constant a = 213 nm. Panel (b): coupling to the second 
mode, a = 294 nm. 



the continuum of radiative PhC modes at the same en- 
ergy. Physically, quasi-guided modes are lossy due to 
out-of-plane diffraction, thus acquiring a finite radiative 
lincwidth 27kn a s discussed in Subscc. IIII Al and calcu- 
lated by the GME method. 

We consider a single QW grown exactly in the mid- 
dle of the GaAs layer, which maximizes coupling to the 
electric field for fundamental slab modes with even par- 
ity. The oscillator strength per unit area will be assumed 
to be f/S = 4.2 x 10 12 cm" 2 , typical of the HH exci- 
ton in an In^Gai-^As QW with thickness £qw = 8 nm 
and low In content The exciton energy is taken to 
be E cx = 1.485 eV and the total exciton mass M ex = 
0.18mo^ For the exciton linewidth we assume a value 
2-f cx = 1 meV. Very high-quality QW structures can 
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Figure 4: (color online) Close-up of Fig. HHb) with exciton 
center-of-mass quantized levels explicitly shown. 
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be currently grown to achieve lincwidths smaller than 
0.5 meV. Here we stay on the conservative side, taking 
into account possible broadening induced by the presence 
of interface defects at the hole boundaries. These consid- 
erations are also supported by preliminary experimental 
data on square lattice PhC slabs with embedded QWs, 
which show no sizeable increase of the exciton linewidth 
as compared to an unpatterned area of the same sample 
(for an initial linewidth on the order of 2 meV)>S With 
the parameters given here, the exciton-photon coupling 
matrix element calculated from Eq. ([50]) is on the order 
of a few meV, with slight variations depending on the 
photonic mode of interest, and the occurrence of strong 
or weak coupling for quasi-guided modes is governed by 
the photonic mode linewidth rather than by the exciton 
linewidth i^i 



A. Guided Polaritons 

Due to the scalability of Maxwell equations, the PhC 
lattice can be engineered in order to have resonance 
between a given QW exciton energy and the desired 
photonic mode at any specific point in the BZ. When 
the exciton-photon resonance occurs below the light line 
(and the exciton linewidth is sufficiently small), the sys- 
tem is always in the strong coupling regime. The new 
quasi-particles that describe the eigenstates of the sys- 
tem are guided PhC polaritons, i.e., their photonic part 
is trapped within the high-index slab through total in- 
ternal reflection. 

In Fig. 05a) and (b) we consider the interaction of a 
QW exciton with photonic modes 1 and 2 from the band 
dispersion shown in Fig. [DJb) , respectively. The photonic 
lattice is engineered to have the resonance condition close 
to the BZ edge along the TX direction. In one case, the 
lattice constant is set to the value a = 213 nm, and in 
the other we assume a = 294 nm, so that the resonance 
at E = 1.485 eV is for k x ~ 0.957r/a [see also arrows in 
Fig-IHb)]. As it can be seen from Fig. [3J the dispersion 
of bare exciton and photon modes is strongly modified in 
both cases, giving rise to sizable anticrossings. It should 
be noted that in Fig. [3j as well as in the following figures 
of Sees. IIVBI and HVCl the uncoupled exciton center- 
of-mass levels are not explicitly shown. We rather pre- 
fer to show the bare exciton and photon dispersions as 
compared to the strongly coupled polaritonic dispersion. 
The calculated vacuum Rabi splitting is hfl^ — 6 meV 
for the first mode and MI-r = 5.5 meV for the second 
one, respectively. The exciton-photon coupling is depen- 
dent on the specific band of interest, due to the different 
spatial profile of the corresponding electric field and thus 
to the modified overlap with the exciton center-of-mass 
wavefunctions. In any case, we point out that such val- 
ues obtained with a single quantum well are compara- 
ble to those commonly achieved for MC polaritons with 
six QWs. 4 The physical reason is the increased exciton- 
photon coupling of Eq. (|30p , due to better confinement 



in the vertical direction of a high-index dielectric slab 
compared to a MC with low-index contrast distributed 
Bragg reflectors. 

A close up of the dispersion diagram showing the 
center-of-mass quantized levels (with M max = 10) is 
given in Fig. QJ for parameters as in Fig. C3Jb). It can 
be noticed that the quantization energy is < 0.1 meV in 
this specific case. These levels represent the effect of spa- 
tial dispersion, well known for bulk exciton-polaritons, 43 
in the presence of center-of-mass quantization. Only the 
first few exciton levels are strongly interacting with the 
photonic mode of interest, due to symmetry of the corre- 
sponding envelope function and electric field. We point 
out that such effects cannot be captured by theoretical 
treatments based on a semiclassical solution of Maxwell's 
equations ; 23 ^ 4 in which the effect of spatial dispersion is 
neglected. Besides giving directly the dispersion of the 
mixed exciton-photon modes, the quantum theory devel- 
oped here is also the starting point for studying polariton 
interactions and nonlinear processes, like in the case of 
bulk and microcavity systems. 



B. Radiative Polaritons 

When a QW exciton is resonant with a quasi-guided 
photonic mode, weak or strong coupling regimes may oc- 
cur depending on the specific situation. We show in Fig. [5] 
the case of a PhC slab of lattice constant a = 430 nm, 
in which the exciton is resonant with different photonic 
modes within the first BZ [namely modes labelled with 
indices 4 and 5 in Fig. [2jb)]. In Fig. [5ja), the bare pho- 
tonic mode dispersion around E c ^ = 1.485 eV is shown. 
It is interesting to notice the existence of a photonic band 
minimum for mode 5 at the T point: such feature leads 
to a quasi-particle dispersion similar to MC polaritons, 
as it will be discussed. Actually, the resonance condi- 
tion occurs simultaneously with different modes along 
TM and TX. In Fig. [5jb), the dispersion of the exciton- 
photon coupled modes is shown in a restricted energy 
range around the exciton resonance. Notice that there 
are five resonant points between QW excitons and PhC 
slab modes in the dispersion diagram, leading to a vari- 
ety of situations for the coupled modes. Along TM, clear 
anticrossings can be seen with photonic modes 4 and 5, 
which are fingerprints of the strong coupling regime. In 
this case, the intrinsic radiative linewidth of bare pho- 
tonic states is lower than the exciton-photon coupling 
energy. As the QW exciton is resonant with mode 4 for 
two different wave vectors along TM, we observe two an- 
ticrossings above the light line in the middle of the BZ. 
This peculiar effect is due to the light dispersion engi- 
neering allowed in PhC structures. Along TX, resonance 
with mode 5 gives strong coupling, while resonance with 
mode 4 at larger in-plane wave vector gives a crossing 
of the bare excitonic and photonic dispersions, meaning 
that the system is in weak coupling. We will analyze 
separately these distinct regimes in the following. 
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Figure 5: (color online) (a) Photonic mode dispersion around the QW exciton energy, E cx = 1.485 eV, for structure parameters: 
r/a — 0.34, d/a = 0.3, a = 430 nm. (b) Close up for the 10 meV energy range of interest, showing the solution for the coupled 
exciton-photon system (full line) together with the bare exciton and photon dispersions (dashed). 



We show in Fig. [5] the complex dispersion of the eigen- 
modes in the region of weak coupling, as compared to 
the bare QW exciton and photon energies. As seen in 
Fig. (Ha), photonic mode 4 is in weak coupling with all 
quantized exciton states (only one is shown here, for 
clarity) close to the X point. A slight modification of 
the exciton dispersion occurs close to resonance for the 
ccnter-of-mass quantized mode having the same symme- 
try as the resonant photonic mode. The photonic mode 
dispersion is not perturbed at all by the presence of the 
QW. Remarkably, the plot of the intrinsic imaginary part 
calculated for the bare quasi-guided photonic mode in 
Fig. [DJb) , evidences a linewidth close to 27 p h ~ 20 meV, 
which explains the weak coupling regime for this particu- 
lar case. Furthermore, the imaginary part of the photonic 
mode energy extracted from the solution of the Hopfield 
matrix shows no modification. On the contrary, and quite 
interestingly, the imaginary part of the QW exciton en- 
ergy, initially set to 7 0X = 0.5 meV for the unperturbed 
state, is sensitively increased around the resonance. This 
indicates the occurrence of an enhancement of sponta- 
neous emission rate or Purcell effect, which may be ob- 
served by time-resolved experiments. 

In Fig. [7| we analyze the strong coupling regime above 
the light line, giving rise to radiative PhC polariton 
states. Mode 5 has a vanishing intrinsic linewidth at 
r [see imaginary part of the bare photonic mode in 
Fig. [7£b), dashed line], leading to a vacuum Rabi split- 
ting. The real parts of the mode energies are shown 
around k = along the TM and TX symmetry direc- 
tions. Looking at the bare QW exciton and mode 5 dis- 
persions, we notice that the resonance condition is not 
exactly at T, but at small wave vectors (|k| ~ 0.027r/a). 
The exact resonant wave vector can be also inferred from 
Fig. Efb) , corresponding to the upper and lower polari- 
ton imaginary parts being coincident and equal to the 
average of bare exciton and photon values. As the bare 



photonic linewidth goes to zero at T, the polariton imag- 
inary part is about half of the bare QW exciton one at 
k = 0. This linewidth averaging effect is well known for 
MC polaritons^ and is a fingerprint of the occurrence 
of a vacuum Rabi splitting. 

C. Discussion and possible experimental 
verification 

PhC slabs allow for a sensitive flexibility and tuning 
capability. In particular, a lithographic tuning (e.g., vari- 
ation of the lattice constant or the air hole radius) on 
different devices patterned on the same chip is commonly 
used to ensure resonance of a desired mode with an ac- 
tive medium. Furthermore, techniques such as digital 
etching^ have been recently employed to achieve a fine 
tuning of photonic mode resonances. By using these post- 
processing techniques, it would be possible to tunc the 
photonic mode of interest with respect to the QW exci- 
ton resonance. In Fig. [5] we show a simulation of such a 
procedure, in which the calculated detunings of the up- 
per and lower polariton branches from the bare exciton 
resonance at k = are reported as a function of hole 
radius. The latter is slightly increased from r/a = 0.336 
(r ~ 144.5 nm) to r/a = 0.344 (r ~ 147.9 nm). The typ- 
ical anticrossing occurs exactly on resonance in the real 
part of energy, and the corresponding imaginary parts 
exchange each other from purely excitonic to purely pho- 
tonic and viceversa. Concerning the significance of Fig. [8] 
with respect to state-of-the art technological capabilities, 
wet chemical digital etching allows for deposition of a 
few monolayers of surface oxide that can be selectively 
removed, thus yielding a wavelength tuning of ~ 1 nm 
per etching cycled On the other hand, surface AFM ox- 
idation of the PhC membrane can lead, in principle, to 
almost continuous tuning of the mode wavelength^ 
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Figure 6: (color online) (a) Real and (b) imaginary parts of 
complex eigenenergies in the weak exciton-photon coupling 
regime close to the BZ edge along the TX direction. The bare 
exciton (dot-dashed) and photon (dashed) complex energies 
are also shown (the real part of photon dispersion nearly coin- 
cides with the coupled system solution). Structure parameters 
are as in Fig. [5] 
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Figure 7: (color online) (a) Real and (b) imaginary parts of 
complex eigenenergies in the strong exciton-photon coupling 
regime close to normal incidence (F-point of the BZ) along the 
two main symmetry directions; the dispersion of uncoupled 
modes is also shown. Structure parameters are as in Fig. [5] 



Radiative PhC polaritons can be probed by angle- 
resolved reflectance from the sample surface, as first 
done on ID PhCs filled with organic molecules^ The 
same kind of experiment could be performed with 
semiconductor-based systems discussed in this work. In- 
deed, the present quantum-mechanical treatment of in- 
teracting photon and exciton states has been compared 
with semiclassical calculations of the surface reflectance, 
showing a very good agreement for the splitting in strong 
coupling regimei 40 ' 41 ! 42 For guided polaritons, on the 
other hand, coupling to an external propagating beam 
is prohibited due to the evanescent character of the elec- 
tromagnetic field in the claddings. In this case it would 
be possible to perform angle-resolved attenuated total 
reflectance by using a high-index prism. This technique 
was applied in Refs. [6511661 to probe the dispersion of pure 
photonic states in line-defect PhC waveguides. 

We now focus on radiative properties close to k = 0. 
As shown in Fig. radiative polaritons can form with 
an energy minimum at the T point, due to the peculiar 
quasi-guided photonic mode dispersion in a square lat- 



tice PhC membrane. The new eigenmodes of the system 
can be excited by coherent or incoherent pumping, and 
probed by emission or reflection at a fixed angle 9, which 
is defined in a vertical plane containing one of the high 
symmetry directions of the PhC lattice. A schematic 
picture of such an experimental configuration is shown 
in Fig. HJa). The angle 9 such that k = (cv/c) sinf?, 
due to conservation of in-plane momentum. We con- 
sider here the angular dispersion along the TX direc- 
tion for the case of zero-detuned (at normal incidence) 
bare exciton and photon modes. The calculated angu- 
lar dispersion is shown in Fig. [9jb) , in which we plot 
the energies of the upper and lower polariton branches 
detuned from the bare QW exciton energy (dot-dashed 
line). It is interesting to notice the similarities between 
the angular dispersion of Fig. [9jb) and the usual dia- 
gram used to illustrate the formation of a polariton trap 
in the energy minimum at k = 0.— The peculiar prop- 
erties of such an angular dispersion have been used in 
the last few years to achieve a number of outstanding 
results requiring non-linear parametric processes of MC 
polaritons . 14 ' 15 ! 16 In particular, coherent population of 
the lowest energy polariton state at k ~ can be realized 
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Figure 8: (color online) (a) Real and (b) imaginary parts of 
upper and lower polariton modes at k — (r-point) as a 

function of hole radius. Structure parameters are as in Fig. [5] -3 -| 1 1 1 
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by polariton-polariton scattering from a precise point in 
the dispersion, which allows for simultaneous energy and 
momentum conservation of the scattered quasi-particles. 
Similar experiments could be realized also with radiative 
PhC polaritons and a parametric process is schematically 
illustrated in the lower polariton branch of Fig. [Hb). 
An extension of the present theory to include nonlin- 
ear terms (such as exciton-exciton scattering)^ in the 
hamiltonian represents a natural extension of this work. 
Moreover, the recent observation of the long-sought Bosc- 
Einstein condensation of MC polariton a 67 ' 68 by use of II- 
VI materials^ has enhanced interest in low-dimensional 
polariton physics. The analogous effect in III-V materi- 
als, like InGaAs/GaAs, has not been observed at time 
of writing. A new route has been suggested to this 
endpS requiring efficient 0D confinement of MC polari- 
tons through the confinement of their photonic part^ 
This is strong motivation for experimental as well as theo- 
retical study of PhC cavity polaritons, in which ultra-high 
quality-factor and small mode volumes can be achieved 
and enhance the radiation-matter coupling in an unprece- 
dented way. 



Figure 9: (color online) (a) Schematic view of angular emis- 
sion from a PhC slab in a plane along a symmetry direction. 
The emission angle, 6 — arcsin(fcc/tj), is defined, (b) Cal- 
culated polariton angular dispersion in the FX direction, at 
zero bare exciton-photon detuning (see Fig. [S| ; a schematic 
polariton-polariton scattering process in the lower branch sat- 
isfying energy-momentum conservation is also shown. 



V. CONCLUSIONS 

In conclusion, we have presented a quantum theory 
to describe radiation-matter coupling for quantum wells 
embedded in a high-index photonic crystal slab with 
a generic pattern. The numerical solution of classical 
Maxwell equations, on which the present theory relies, 
has been previously reported within a guided-mode ex- 
pansion approach^ which yields not only the photonic 
mode energies but also the crucial photonic linewidth 
parameter. After obtaining the second-quantized total 
hamiltonian of the system in the linear regime of low 
excitation density, we have diagonalizcd it through a 
generalized Hopfield method. Thus, the complex eigen- 
modes of the exciton-photon coupled system have been 
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obtained. Both weak and strong coupling regimes are 
treated within the present theory, allowing for an effi- 
cient design of photonic lattices for specific purposes. 

We have shown results on the square lattice PhC slab 
demonstrating the formation of guided and radiative pho- 
tonic crystal polaritons. The latter are formed when 
the exciton-photon coupling exceeds the intrinsic photon 
linewidth. In strong coupling, the vacuum Rabi splitting 
at resonance is larger than for the corresponding micro- 
cavity polaritons, due to the increased field confinement 
and better overlap of exciton and photon wave fields in 
a PhC slab. As a consequence, a more robust polari- 
ton effect can be envisioned in such structures, which 
could be useful for nonlinear polariton applications. Re- 
lated to the latter point, we reported on a specific lattice 
design for which the radiative polariton dispersion has 
a minimum around normal incidence. Such dispersion 
closely mimics the one obtained in usual microcavities, 
which acts as a polariton trap and is at the origin of 
current research on nonlinear parametric processes with 
exciton-polaritons. We believe that the present results 
will stimulate further research on non-linear polariton ef- 
fects in semiconductor-based photonic crystals, which in 
turn could connect with polariton quantum optic a 71 ' 72 ! 73 
as an emerging field of research. 



a Bogoljubov transformation on bare exciton and pho- 
ton operators. It was originally used by Hopfield for the 
case of bulk polaritons^ and then generalized to polari- 
tons in planar MC o 30 ' 34 and in micro-pillars ^ Here we 
provide a detailed derivation of the generalized Hopfield 
matrix for the case of PhC slabs, which has to be diag- 
onalized numerically The second-quantized total hamil- 
tonian, Eq. (|32p . and the polariton operator expansion, 
Eq. (|33|) . are substituted in Eq. (|35| . The general com- 
mutation relations for bosonic operators 



#k,k'<!>n,n'<ikn 

— <5k,k"5n,n'a kn 
<5k,k'<5n,n'akn 

— ^k.k'^ra.n'Skn 



(Al) 



Appendix A: GENERALIZED HOPFIELD 
METHOD 

The Hopfield method to diagonalizc the coupled 
exciton-photon hamiltonian is equivalent to performing 



must be satisfied. Thus, after factorizing terms with com- 
mon operators, the following system of linear equations 
in the variables w n , y n , x v , and z v can be derived (we 
drop the fixed subscript k rcsp. — k, for easier notation): 



hui n w n + ^ ' 2D unn > w n > — ^ 2D unn i y n r — S ' iC nv x v — S ' iC m z v — HQ w n 

+ ^ iCm Wn ~ X! iCw Vn + ^ cx) Xv = ^ Xn 

n n 

^ ^ nn' l^n' h(jJ n y n ^ ^ ^Dy nn' Vn' ^ ^ ^C n y X v — ^ ^ iCnv %v — HQ y n 

n' .v n' ,v v v 

- ^2 iG ™ w n + ^2 iCnv Vn ~ £^° x) z v = hVlz n . (A2) 

n n 

i 



The sums in Eq. (|A2j) must be truncated in order to deal with n = 1, . . . , iV max , v = 1, . . . , -A/ max . The explicit 
with finite matrices. If iV max photonic bands and Af max form of the generalized Hopfield matrix is given by 
excitonic levels at fixed k are retained in the expansion 
\, the problem is reduced to an eigenvalue equation 



Mv = Mlv , (A3) 

where the matrix M has dimension 2(A r max + M max ) x 
2(A r max + M max ), and it is a generalization of the 4x4 

matrix derived by Hopfield^ The vector v is simply M = I "^,t^ I • (A5) 

—I— ^ I I — it. — / . t — ^ I I — t I . * ' 

v= (w n ,x v ,y n ,z v ) T , (A4) 
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The single blocks are respectively given by the -ZV max x 
AL™ matrices 
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w + 2D=[fej n <y n , n /+5^2D vnn# ], (A6) 



by the M max x M max diagonal matrices 

E = \E^5 uy \ , 
and by the iVmax x M max matrices 
C = [C n „], 



(A7) 



(A8) 



while denotes the M max x M max matrix with zero en- 
tries. It should be noted that both positive- and negative- 
defined energy values arc present in the Hopficld matrix. 
After the numerical diagonalization, yielding the eigen- 
values http = E p +i~f p , with p = 1, 2, . . . , 2(7V max +M max ), 
only those with positive real part are eventually retained. 
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